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ABSTRACT 

The angular diameters of six oxygen rich Mira-type long-period variables have been measured at var- 
ious near-infrared (NIR) wavelengths using the aperture masking technique in an extensive observing 
program from 1997 Jan to 2004 Sep. These data sets span many pulsation cycles of the observed objects 
and represent the largest study of multi-wavelength, multi-epoch interferometric angular diameter mea- 
surements on Mira stars to date. The calibrated visibility data of o Cet, R Leo, R Cas, W Hya, x Cyg 
and R Hya are fitted using a uniform disk brightness distribution model to facilitate comparison between 
epochs, wavelengths and with existing data and theoretical models. The variation of angular diameter 
as a function of wavelength and time are studied, and cyclic diameter variations are detected for all ob- 
jects in our sample. These variations are believed to stem from time-dependent changes of density and 
temperature (and hence varying molecular opacities) in different layers of these stars. The similarities 
and differences in behaviour between these objects are analyzed and discussed in the context of existing 
theoretical models. Furthermore, we present time-dependent 3.08 /mi angular diameter measurements, 
probing for the first time these zones of probable dust formation, which show unforeseen sizes and are 
consistently out of phase with other NIR layers shown in this study. NIR light-curves were recovered, 
and show the distinctive phase lag in the m axi ma of » 0.2 (compare d to the visual maximum), similar 



to the lag found by, e.g. ISmith et all d2002h and lNadzhip et al.l d2Q0lb . The S-type Mira x Cyg exhibits 



significantly different behaviour compared to the M-type Miras in this study, both in its NIR light-curves 
and its diameter pulsation signature. Our data show that the NIR diameters predicted by current models 
are too small and need to incorporate additional and/or enhanced opacity mechanisms. Also, new tailored 
models are needed to explain the behaviour of the S-type Mira^- Cyg. 

Subject headings: instrumentation: interferometers - techniques: interferometric - stars: late-type - stars: AGB and 
post-AGB - stars: fundamental parameters - stars: individual: Mira stars: individual: R Hya stars: individual: Chi 
Cyg stars: individual: W Hya stars: individual: R Leo stars: individual: R Cas 
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length and pulsation phase (e.g., lHaniffetal.il 19951: 
van Belle et ail 1 1 9961: Pen-in et all fl999l: I Young et al I 
IQOllHofmann et al. 120021: iThompson et all l2002t 
Woodruff et al. 2004). Interferometric studies of Mira 



stars conducted over multiple epochs and spanning a 
range of wavelengths have helped to address funda- 
mental que stions such as the pulsation mode of thes e 
stars (e.g., I Woodruff et al. l2004HFedele et al l 12005b . 



the m olecular and dust abu n dances in the atmospher e 
(e.g. Ilreland et all l2004albl: llreland & Scholj 120061) . 
the characteristics of t he circumstellar environment 
(e.g. lDanchi et a l. 1994) and photosphe ric/circumstellar 
asymmetries (e.g. iRagland et al.ll2006l) . 

Interferometry, together with a host of other obser- 
vational techniques, has gradually been making ad- 
vances into our understanding of the basic physics of 
these stars. For example, there is now consensus that 
Miras pulsate in the fundamental mode, bas ed on pho- 
tomet ry, e.g., MACHO obser vations (see | Wood et al.l 
Il999h ." and spectroscopy (see IScholz & W ood 2000) 
in addition to interferometry. However, the struc- 
ture and composition of the stellar atmosphere still 
shows surpris i ng re sults when investigated in detail. 
Ragland et al. (2006) have, e.g., detected H band 



asymmetric brightness distributions in about 29% of 
their sample of nearby asymptotic giant branch stars, 
substantiating another level of complexity in the struc- 
ture of this class of objects. 

In this paper we present the most comprehensive in- 
terferometric study of Miras to date, encompassing 6 
nearby objects, observed at up to 19 different phases in 
4 filters. With this homogeneous data-set, we are able 
to examine phase-dependent variations in Miras' at- 
mospheric structure and are in a position to investigate 
long term effects that span several pulsation cycles. By 
observing in the NIR (especially within the J, H and K 
bandpasses), we can sample molecular strata near the 
continuum-forming layers that are often clos e to the 
positi on of the Rosseland layer (t Ross = 1, cf. Schold 
2003 and references therein). 



Theoretical models describe the molecular layer- 
ing in the atmosphere and its variation with time (e.g. 
Bessell et all 1 9961: [kofmann et al.ll 19981: llreland et al 



200413)7 and make predictions of observables such 
as light-curves, diameters and intensity distributions. 
With our extensive database we can challenge existing 
models and motivate the development of a more com- 
plete physical picture of Mira variables. Ultimately, 



the aim of interferometric observations like these is 
to calibrate the theoretical models so that fundamen- 
tal parameters of Miras can be derived from simple 
observables, enabling studies of stellar populations be- 
yond the solar neighbourhood. 

2. Observations and Data Reduction 
2.1. Aperture-masking Observations 

Our sample of 6 Miras contains only nearby ob- 
jects, all of which have revised Hipparcos distances 
(see Table [U. The objects were chosen for their large 
angular diameters and NIR brightness. Observations 
were performed with the 10 m Keck I telescope at 
a range of NIR wavelengths (see Table |2]i using the 
Near Infrared Camera (NIRC). The telescope pupil 
was converted into a sparse interferometric array by 
placing aperture masks in the beam in front of the in- 
frared secondary mirror, allowing the recovery of the 
Fourier amplitudes and closure phases for baselines 
up to 9.8 m. For a detailed discussion of mask design, 
observing methodolog y, scientific rationale and imple- 
mentation, we refer to lTuthill et al. ( 2000bl) . 



For this work we used non-redundant masks with 
15 or 21 holes configured to deliver near-o ptimal sam- 
pling of the Fourier plane (IGolayi 1 1 97 lh . Data-sets 
consisting of 100 140 msec exposures were taken, al- 
ternating between the target of interest and nearby cal- 
ibrator stars. The latter were chosen to have well char- 
acterized, smaller apparent sizes (see Table |3). The 
data were recorded at 19 different epochs spanning 
more than 7 years, delivering good coverage through 
the pulsation cycles of the objects observed. Tables |4] 
to|9]list these observations. 

2.2. Extraction of Visibilities 

The procedures for extracting the visibility ampli- 
tudes, as well as en gineering and pe rfor mance details, 
are do cumented in Mo nniei ( 1999b and Tuthill et al 



(2000b), while recent scientific app lications o f the 
data pipeline can be found in lMonnier et al .1 (12002b and 
Tuthill et al.l (12002b . In principle, the pupil geometry 
of the telescope mimics the operation of a separate- 
element interferometer array, and the data collection 
and analysis are similar to standard methods for inter- 
ferometry experiments such as speckle imaging. The 
short-exposure images are dark-subtracted, flat-fielded 
and cleaned of pattern noise. Power spectra are then 
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Table 1 
Observed objects 



Name Period 


Spectral 


Distance 


Distance 


[days] 


Type Range 


[pc] 


reference 



o Cet 


332 


M5-9e 


107 ±6 


1 


RLeo 


312 


M6-9.5e 


82 + 5 


1 


RCas 


430 


M6-10e 


100 + 5 


2 


WHya 


385 


M7.5-9ep 


78 + 3 


1 


*Cyg 


408 


S6-S9(MSe) 


149+11 


1 


RHya 


380 


M6-9eS(Tc) 


118 + 7 


1 



References. — Object Period from the American Associa- 
tion of Variable Star Observers (AAVSO) visual light curves 
(A.A. Henden et al. 2 006, private commun ication), M Spec- 
tral Type Range from Sloan & Price ( 1998) and^ Cyg Spec- 
tral Tyrje_fromjKeenan & Boeshaa rl dl980h. Distan ces from 
(l)lKnapp et all d2003l) or (2)|Pourbaix et alj d2002l) 



Table 2 
Table of filters 



Filter Keck/NIRC Centre Wavelength Bandwidth 
Filter Name [/vm] [//m] 



zl.08 
J 1.24 
H1.65 

K2.26 

L3.08 
L3.31 



Hel 
Oil 
Fell 
Kcont 
H221 
PAHcs 
PAH 



1.083 
1.236 
1.647 
2.260 
2.261 
3.082 
3.310 



0.014 

0.011 

0.018 

0.050 

0.0239 

0.101 

0.063 



3 



Table 3 

Calibrator stars with estimated diameters 



Calibrator Spectral 


Adopted UD Reference 


Type 


Angular Diameter [mas] 



aCet 


M1.5III 


11.7 + 0.6 


1 


a Ari 


K2III 


5.9 + 0.6 


1 


a Cas 


Kom 


5.3 + 0.1 


2 


/?Cet 


Kom 


6.7 + 0.5 


4 


a Hya 


K3II-III 


9.1+0.1 


3 


n Leo 


M2III 


4.6 + 0.3 


5 


a Lyn 


K7III 


7.2 ± 0.6 


1 


2Cen 


M4.5IH 


14.7 


6 


7 Sge 


MOIII 


6.0 + 0.6 


1 


<5Sgr 


K3III 


6.9 + 0.9 


7 


4 Cas 


Ml III 


4.3 


6 


fCyg 


KOIII 


4.3 + 0.4 


2 


(JOph 


M0.5III 


10.1+0.5 


8 


Vega 


AOV 


3.1+0.1 


9 


£Cyg 


K4.5I 


7.5 + 0.6 


10 


y Hya 


G8III 


3.0 + 0.2 


11 


7r Hya 


K2III 


3.9 + 0.2 


4 



Ref erenc es. — (1) Dycketal. (1998); (2) iMozurkewich et al 
Jl99ll): (3) iMozurkewich et all J2003 b; (4) from CHARM cata- 
log dRichichi& Percheronll2002h: (51 iRidewav et all dl979b: (6) 
iDumm & Schildl dl998b: (7)lMonnier et al.l ( 12004: 78)|Perrin et all 
dl998h : (9) lHanburv Brown et all dl974 : ( lOt iDvck et al.l dl996l) : 
(11) from CHARM2 catalog dRichichi et al.ll2005h 
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Table 4 
Observations of o Cet 



Date 


JD 




UDji.24 


UDhi.65 


UD K 2.26 


UD L 3.08 


Calibrators 




-2450000 




[mas] 


[mas] 


[mas] 


[mas] 




1997Decl6 


800 


0.94 


22.2±1.0 


27.5 + 2 2 2 4 


31.0+1.9 


58.8+1.5 


a Cet 


1998Sep29 


1056 


1.71 


29.1 ±2.2 


33.3+1.6 


37.0+1.2 


57.7+1.1 


a Cet 


1999Jan05 


1184 


2.10 


25.6±1.3 


27.6 + 2 2 f 


31.8+1.8 


61.4+1.3 


a Cet 














61.5+1.2 
















61.5+1.3 
















61.5+1.2 




1999M29 


1390 


2.72 


29.4±1.9 


32.9+1.7 


36.9+1.2 


53.2+1.4 


a Cet 


2000Jan25 


1570 


3.26 




30.1+2.6 


33.1 + 1.3 


58.3+1.4 


a Cet 


2000Jun23 


1720 


3.71 


29.9±1.9 


38.1 + 1.8 




56.9+1.3 


a Ari 














58.2+1.0 




2001M29 


2121 


4.92 




30.1 + 1.9 


33.8+2.7 


58.7+5.0 


a Cas , a Cet 


2002Jul23 


2479 


5.98 


24.8±2.1 


26.6+0.5 


34.9+1.5 


61.5+2.3 


a- Cet 


2004Sepl4 


3262 


8.30 


27.1±1.0 


29.3+1.6 


32.3+1.2 


59.0+1.4 


a Cet , /3 Cet 












31.8+0.9 


57.5+1.2 





57.4+1.1 
58.2+1.9 
57.8+1.1 
57.9+1.4 
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Table 5 
Observations of R Leo 



Date 


JD 


O 


UDji.24 


UDhi.65 


UD K 2.26 


UDL3.08 


Calibrators 




-2450000 




[mas] 


[mas] 


[mas] 


[mas] 




1997Jan29 


478 


0.05 


31.0+1.1 


29.6+1.9 


30.3+2.6 


52.1 + 1.5 


a Hya , n Leo 








32.6±1.5 




32 ' 9 -2 2 2 


52.8+1.4 
















55.6+1.7 
















54.6+1.4 
















53.8+1.2 
















54.5+1.5 




1997Decl6 


800 


1.14 


30.0+1.1 




31.4+1.7 


47.8+1.8 


a Hya 








29.9±1.2 










1998Aprl4 


918 


1.54 




32.9+1.4 


32.3+1.8 


49.2+1.6 


a Lyn 


1998Jun04 


970 


1.71 


33.1±2.0 


29.7+2.0 


32.6+1.9 


50.7+1.6 


a Lyn 


1999Jan05 


1213 


2.40 






33.9+0.9 




n Leo 


1999Feb04 


1213 


2.49 


31.2+1.8 


33.2+2.0 


35.4+ 3 9 3 


50.9+1.9 


a Hya 


1999Apr25 


1295 


2.75 


29.6+1.4 


29.1+0.8 


34.6+1.3 


53.3+1.5 


a Lyn, n Leo 












34.1 + 1.4 






2000Jan25 


1570 


3.64 




34.6+1.6 


33.1 + 1.7 


52.6+0.7 


a Lyn 














52.1+0.8 




2000Jun23 


1720 


4.12 




29.3+2.1 


31.2+1.9 


55.4+1.8 


7rLeo 


2001Junll 


2073 


5.24 


31.5+1.3 


33.5 + 3 4 8 " 


30.7+ 4 8 5 


56.5+2.0 


n Leo , a Lyn, 2 Cen 


2003Mayl2 


2772 


7.48 




36.2+1.1 


37.r 4 5 4 2 




7rLeo 


2004May28 


3154 


8.71 




30.1 + , 2 8 ° 


32.1 + 1.7 


56.4+1.4 


n Leo 


2005May26 


3516 


9.88 




28.9+2.0 


31.5+1.5 




7rLeo 
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Table 6 
Observations of R Cas 



Date 


JD 




UDji.24 


UDhi.65 


UDk2.26 


UDl3.08 


Calibrators 




-2450000 




[mas] 


[mas] 


[mas] 


[mas] 




1997Dec08 


802 


0.67 






30.0+1.5 




a Cas 


1998Jun04 


970 


1.06 


25.0±1.2 


23 .2+ 3 , 1 


24.2+1.6 


44.8+2.2 


y Sge, a Cas 


1998Sep29 


1056 


1.27 


24.5+1.5 


27.0: 2 2 2 4 




41.1+3.0 


a Cas 


1999Jan05 


1184 


1.57 


24.4+1.8 


23.9_; 2 4 8 


28.2+1.2 


44.1 + 1.8 


a Cas 


1999M29 


1390 


2.04 


24.1 + 1.8 




28.5+1.4 


48.3+1.5 


a Cas 


2000Jan25 


1570 


2.46 




27.7+ 2 2 


29.7+1.2 




a Cas 


2000Jun23 


1720 


2.81 


24.9+1.4 


26.1-f 


28.9+1.9 


46.2+1.6 


a Cas 


2001Junll 


2073 


3.63 




26.5+ 2 2 4 


27.4+ 4 54 




a Cas 












27.6£ 4 
















27.7_ + 4 5 3 2 






2001Jul29 


2121 


3.74 




25 5 +2J 

* -2.4 


26.1+2.1 




a Cas 


2002Jul23 


2479 


4.58 




29.5+0.9 


29.6+1.5 


42.0+2.0 


a Cas 










27.4+0.9 








2003Mayl2 


2772 


5.26 




23.3+2.3 


22.1 + 7 143 




5Sgr 










24.4+1.4 


29.0- 






2004May28 


3154 


6.14 




23.2 + 2 2 5 5 


24.6- 2 


48.9+1.5 


4 Cas 


2004Sep2 


3262 


6.39 


27.6+0.9 


26.8- 


25.2+1.5 


42.2+1.3 


P Cet , a Cas ,e Cyg 








26.6+1.1 


26.5+1.8 


24.5+1.3 


40.5+2.1 












27.9+2.0 
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Table 7 
Observations of W Hya 



Date 


JD 


<5 


UDj i.24 


UDhi.65 


UDk2.26 


UDl3.08 


Calibrators 




-2450000 




[mas] 


[mas] 


[mas] 


[mas] 




1997Jan29 


478 


0.53 


41 9 +32 

-3.2 


38.9+1.8 


42.4+1.8 


62.3+1.4 


£Oph 








42.0+ 4 ., 1 


38.3+1.0 


42.5+2.2 


61.7+1.1 














41.9+1.8 


60.3+1.0 














39.4+2.0 


60.6+1.5 














39.7+2.7 
















41.4+2.2 






1998Aprl4 


918 


1.76 


32.8+2.4 


34.2+1.7 


37.4+1.8 


62.3+1.2 


a Lyn, 2 Cen 










34.6+2.0 


39.8+1.4 
















36.7+1.2 






1998Jun04 


970 


1.90 


32.0+2.3 


33.6+1.8 


36.4+1.1 


62.7+1.0 


2Cen 


1999Jan05 


1184 


1.50 






43.0+1.4 




2 Cen 












41.8 +1.2 
















/IO A _i_1 7 
















42.6 +1.8 






1999Feb05 


1213 


2.58 


39.3+3.0 


41.6+1.2 


42.0+ 2 , 7 


63.2+1.2 


2 Cen 












43.8t, 2 2 4 






1999Apr25 


1295 


2.79 


35.2+3.5 


38.1+2.4 


41.0+1.0 




2 Cen 










36.7+2.2 


36.6+1.6 














36.3+1.7 
















36.8+1.8 








1999M29 


1390 


3.04 




34.8+1.6 


36.7+0.9 




2 Cen 


2000Jan25 


1570 


3.53 


47.1+2.7 


43.2+1.5 




64.5+1.1 


2 Cen 












45.9+0.8 






2000Jun23 


1720 


3.94 




35.0+1.8 




68.5+1.6 


2 Cen 


2001Junll 


2073 


4.88 


34.5+3.7 


36.0+1.8 


41 7 +3J 


65.3+1.4 


a Lyn, 2 Cen 












39.4^' 






2003Mayl2 


2772 


6.66 




40.6+1.3 


44.2+ 3 3 8 




2 Cen 


2004May28 


3154 


7.56 




39.1 + 1.5 


41.5+1.4 


65.7+1.2 


2 Cen 












41.4+1.2 






2005May26 


3516 


8.55 




38.4+1.9 


43.8+1.4 




2 Cen 
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Table 8 
Observations of \ Cyg 



Date 


JD 


O 


UDji.24 


UDhi.65 


UD K 2.26 


UD L 3.08 


Calibrators 




-2450000 




[mas] 


[mas] 


[mas] 


[mas] 




1998Jun04 


970 


0.59 


21.3+2.8 


19.5 + 3 3 f 


23.0+2.1 


45.6+3.0 


y Sge 


1998Sep29 


1056 


0.87 




2i.2_ + , 3 9 4 


26.1+2.7 


46.0+4.0 


y Sge 


1999Jul29 


1390 


1.62 


21.6+1.6 


22.4+ 3 7 2 


28.2+1.5 


42.4+3.0 


y Sge, Vega 














43.7+2.8 




2000Jun23 


1720 


2.48 




23.1^ 


25.1+2.2 


39.5+2.4 


7 Sge, £ Cyg 












24.9+ 2 f 


41.0+2.7 




2001Junll 


2073 


3.36 


25.4+1.9 


26.4+ 28 


30.4; 4 ; 7 


45.5+2.5 


y Sge 


2002Jul23 


2479 


4.36 


28.5+2.7 


27.5+1.0 


29.6±1.6 


44.3+2.5 


y Sge 


2003Mayl2 


2772 


5.07 




28.1+2.0 


49 1 +4-4 


46.5+2.6 


eCyg 


2004May28 


3154 


6.01 




24.7+ 3 , 


33.8+3.5 


53.1+4.0 


Vega 


2004Sep23 


3262 


6.31 


24.8+1.2 


26.7+2.0 


33.0+2.3 


50.5+3.9 


y Sge 



Table 9 
Observations of R Hya 



Date 


JD 


O 


UDji.24 


UDhi.65 


UDk2.26 


UD L 3.08 


Calibrators 




-2450000 




[mas] 


[mas] 


[mas] 


[mas] 




1997Jan29 


478 


0.70 


26.2+1.9 




26.0: 28 


40.0+2.1 


^Oph 


1999Jan05 


1184 


2.56 






24.3+2.3 
27.7+0.9 
27.4+1.0 
24.0+1.9 




y Hya 


1999Feb05 


1213 


2.64 




30 +2 2 

-'"•"-2.0 






n Hya 


1999Apr25 


1295 


2.86 




24.1+2.9 


28.0+2.2 




2Cen 


2000Jan25 


1570 


3.58 






31.0+1.3 




2Cen 


2003May28 


2772 


6.74 






27.2^ 




2Cen 
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computed frame by frame as the squared modulus of 
the Fourier transform. Stellar fringes appear as dis- 
crete peaks in such power spectra, with the origin 
occupied by a peak whose height is proportional to 
the squared flux in the frame. Squared visibilities are 
found by dividing the power at the spatial frequency 
of the fringes by that at the origin and then normaliz- 
ing with the corresponding signal from the calibrator 
spectrum. The uncertainty associated with the squared 
visibilities is derived from the scatter in each ensemble 
of 100 exposures. 



2.3. Seeing Correction 

For all baselines and for all target and calibrator 
stars considered in this paper, the dominant noise 
sources were seeing and windshake (wind-induced 
telescope wobble). When the atmospheric conditions 
vary between observing the source and its calibrator, 
the overall ratio changes between the fringe power 
and the total flux on the detector. As clarified by 



Monnier et al.l (2004), for aperture-masking data this 
change is nearly constant as a function of baseline for 
baselines longer than the coherence length ( w 0.5 m 
at K2.26 band). This means the visibility function 
will approach a non-unity value at short baselines. We 
took the visibility at the origin to be a free parameter 
in our studies, an assumption that does not affect our 
fitting procedures as long as there is no significant flux 
coming from over-resolved structures (> 0".5). See- 
ing and windshake also lowered the mean visibility 
V and increased its variance on each baseline. How- 
ever, the quantity was independent of the stellar 
brightness and the degree to which it was resolved 
by the baseline in question, only depending on seeing 
and windshake. We partially corrected for the effects 
of seeing and win dshake using the method described 
bv llrelandl d2006l) . which entails empirically fitting a 
function of the form V — Voexp(-fc ^ a ^ ) ) to indi- 
vidual star observations and correcting visibilities for 
each target and calibrator star prior to dividing the tar- 
get with calibrator visibilities. We used a conservative 
value for k of 0.8 (which may have under-corrected for 
the effects of seeing and windshake) and smoothed the 
function ^2 m t h e u — v plane in order to minimise 
err ors in applyin g this correction (for further details 
Ireland! 120060 . We evaluated the effectiveness of 



The apparent angular diameters were closer to the true 
values after the seeing/windshake correction on > 70% 
of the stars examined (see FigureHJ. Only on 3 object 
observations with very bad windshake did the algo- 
rithm cause a marginal worsening of the miscalibration 
(5 mas). 

Figure [2] shows typical data obtained from the 
masking experiment. The targets were seldom fully 
resolved, making image reconstruction only possible 
for those objects with the largest angular diameters 
(e.g. W Hya). These results will be presented in a 
subsequent paper. The two dimensional visibilities 
showed no significant deviation from circular symme- 
try (except for W Hya), possibly due to the lack of high 
spatial frequencies sampled by our <10m baselines. 
Since W Hya's deviation from circular symmetry is 
smaller than the UD angular diameter error, we made 
the assumption that all the stars are spherical and the 
data were azimuthally averaged. 

To keep the representation of the data homoge- 
neous, our quantitative analysis is based on fitting sim- 
ple uniform disk (UD) models to the azimuthally av- 
eraged Fourier data. We concentrate on the four fil- 
ters which provide most of the temporal coverage in 
our observations, namely L3.08, K2.26, H1.65 and 
J 1 .24, considering the other filters only in Section l4~T1 
Here, e.g., L 3.08 indicates that the filter covers a sub- 
interval of the L bandpass of the conventional UB V fil- 
ter system and is centered at 3.08 /im (Table 13. Note 
that the K 2.26 bandpass encompasses two very similar 
filters, both centered at 2.26 fim with slightly different 
bandwidths. 

Calibrator stars nearby in the sky and measured 
interleaved with the target observation are used to 
estimate the system transfer function, a standard 
practise in interferomet ry described in detail by 



see 



this correction by comparing the calibrated visibilities 
of unresolved stars with and without seeing correction, 
which turned out to improve the calibration process. 



iMillan-Gabet et al.l (12005b . 

The visibility amplitude was calibrated by dividing 
the target visibilities by the calibrator star visibilities, 
after first correcting for the estimated sizes of the cali- 
brators (Table |3]l. 

2.4. Uniform Disk diameters 

Although the true stellar intensity profile is not a 
UD, fitting the observed visibilities with this simple 
profile still provides a useful estimate of the apparent 
size of the target. Our longest baseline was just under 
10 m, and so we are resolving low-resolution structure 
in the target star's intensity profile. This makes it dif- 
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Fig. 1. — Left: Calibrated visibilities (without azimuthal averaging) of an unresolved target star (a Ari). Due to the 
effects of variable seeing and/or windshake, the object appears to have an UD angular diameter of 23 mas. Right: The 
same data with the seeing correction applied, showing that these data are now consistent with an unresolved source. 
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o Cet at Phase 0.71 



□ 3.08um Diam = 57.7 +/- 0.4 mas 

A2.26um Diam = 37.0 +/- 1.2 mas 

;> 1.65'um Diam = 33.3 +7- 1.6 mas 

* 1.24/im Diam = 29.1 +/- 2.2 mas 




R Leo at Phase 0.49 



D3.0Bum Diam = 50.9 + /- 1.9 mas - 

. ■ i" Diam = 35.4 +/- 3.3 mas 

C 1.6;";.v.in Dinm = 33.2 +/- 2.0 mas 

* 1.24^.m Diam = 31.2 +/- 1.6 mas - 




2x10" 3xlO b 4x10^ 

Spatial Frequency (rad -1 



W Hya at Phase 0.76 



□ 3.08/j.m Diam = 62.3 +/- 1.2 ____ 

i I i m = 37.4 +/- 1.8 mas - 

■> 1.65;,-in_ Diam = 34.2 + /- 1.7 mas 

* 1.24/xm Diam = 32.8 +/- 2. * 



2xlO u 3xl0 u 4xlO u 

Spatial Frequency (rad"' 
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X Cyg at Phase 0.6S 



□ 3.08um Diam = 43.7 +/- 2.8 mas 

A2.26um Diam = 28.2 +7- 1.5 mas 

O 1.65um Diam = 22.4 +/- 2.9 mas 

* 1.24/im Diam = 21.6 +/- 1.6 mas 




R Hya at Phase 0.70 



2xl0 6 3x10 s 4xl0 6 

Spatial Frequency (rad" r 



□ 3.08um Diam = 43.7 + /- 2.8 mas 
i2.26 ;j .m Diam = 2 8. 2 +7- 1.5 mas 
O 1.65um Diam = 22.4 +/- 2.9 mas 




2xlO e 3xlO fl 4xl0 fl 

Spatial Frequency (rad"') 



Fig. 2. — Azimuthally averaged squared visibility functions for o Cet, R Leo, R Cas, W Hya, x Cyg and R Hya at 
various visual phases Results for the different filters are offset for clarity. The very large apparent diameters in the 
L 3.08 filter can be attributed mainly to the H2O opacity, as first suggested and observed bv lTuthill et all d2000al) . The 
diameters correspond to best-fit UD models. 
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ficult to differentiate between a UD, a fully-darkened 
disk, a Gaussian or a more complex intensity distri- 
bution. We chose UD diameters to allow compari- 
son of findings with existing literature and to avoid 
the difficulties encounter ed when deriving more so- 
phisticated diameters (cf. iHofmann et al ] ll998HScholj 
2003). If the "true" intensity distribution is known, or 
a predicted model distribution is to be compared with 
the data, different radius definitions can be easily con- 
verted from and to a UD radius. 

In various stars, the J 1 .24-band curves show some 
deviation from the UD model at around 3 x 10 6 rad _1 , 
which could indicate the presence of dust causing scat- 
tering in the outer layers, leading to a smaller "true" 
photospheric diameter (see Figure|2]i than derived with 
a simple UD model fit. In addition, due to their larger 
angular diameters (especially in the L3.08 filter), W 
Hya and possibly o Cet were also resolved enough to 
detect further deviations of the intensity distribution 
from the simple UD shape to be compared with model 
predictions. They are also possible targets for imaging 
after recovering the closure-phase information in ad- 
dition to the Fourier amplitudes. These comparisons 
with more complex models and interpretations will be 
addressed in a future paper. 

2.5. Estimation of Seeing Miscalibration 

Given that the seeing correction described in Sec- 
tion [23J leaves some residual error and assuming that 
the dominant error term is the uncertainty introduced 
by variations in the seeing between measurements of 
targets and calibrators, we undertook a separate study 
which assessed the robustness and repeatability of the 
calibration process to seeing-induced miscalibration. 
If the seeing remains constant, then calibrating (i.e., di- 
viding) the visibility function from successive datasets 
taken on the same object should yield the visibility 
function of a point source (i.e V 2 = 1). In the case 
of changing atmospheric conditions, considerable dif- 
ferences in the transfer function can appear between 
consecutive datasets. The error induced by such fluc- 
tuations can be measured from observations of a sin- 
gle object (either target or calibrator) taken through the 
same filter at different times during one night. 
Once again, we chose to fit a UD model to the data, 
a natural candidate to fit the error introduced by the 
change in seeing throughout one observing night. This 
allows us to estimate atmospheric uncertainties by ap- 
plying the same fitting procedure as was used for the 
rest of our data, reducing the risk of biases introduced 



by the fitting method. 

Asymmetric errors were calculated by calibrating 
the target with a UD of the derived diameter. These 
errors pertain to the respective filters, and if no object 
was observed twice in one night in the same band, the 
maximum global error is assumed. 

The data in Figure [2] demonstrate another well 
known problem: that of calibration near the so-called 
seeing spike. Note how the azimuthally averaged vis- 
ibility points at short baselines can deviate consider- 
ably from the otherwise uniform shape of the visibility 
curve. As the calibration of these large fluctuations at 
low spatial frequencies is highly challenging, we fitted 
our models only to the spatial frequencies in Fourier 
space corresponding to baselines of 2m or more. 

2.6. NIR photometry 

In the course of the data reduction, we extracted the 
total received flux for both science objects and cali- 
brators. With this we were able to retrieve contem- 
poraneous photometry data in the the same bands as 
our angular diameter measurements. Thes e measure 



merits complement the photometry by lWhitelock et al. 



(2000), which does not cover our full range of ob- 
servations, and provide a powerful tool to monitor 
NIR light-curve s . The small differences between the 
Whitelock et al.1 d2000l) photometry and our measure- 



ments can be attributed to our use of filters of differing 
bandpass and center wavelength(see Table|2|l. 

3. Results 

3.1. Light-curves 

Figure [3] shows our NIR photometry plotted as a 
function of pu lsation phase, togethe r with the JHK 
photometry of lWhitelock et all d200d) and AAVSO vi- 
sual photometry (A. A. Henden et al. 2006, private 
communication). As an approximation based on the 
appearance of the light curves, we fitted our NIR pho- 
tometry data with simple sinusoidal functions (not 
shown in the figures), with the exception of R Hya, 
where the light curve was too incomplete. Although 
these fits are only approximations, they provide a use- 
ful tool to examine the dates of maxima and minima 
as a function of wavelength. The NIR magnitude vs. 
phase curves are certainly not strictly sinusoidal, but 
depending on the star and the bandpass they can be ap- 
proximated by a sine function. Our photometry shows 
a small lag of up to 0.07 cycles in the maxima of 
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Fig. 3. — All Figures: Our NIR photometry folded with the visual pulsation period, shown with Wh itelock et al. (2000) 
JHK photometry and visual light-curves. The J 1.24,H 1.65,K2.26 and L 3.08 fluxes extracte d from Nirc masking ar e 
represented by red asterisks, blue diamonds, green triangles and black squares respectively. Wh itelock et al. hood J 



(red crosses), H (green crosses) and K photometry data (red crosses) are shown where available. Visual data kindly 
supplied by the AAVSO. 
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light curves from the J 1.24 to the H 1.65 and from the 
H 1.65 to the K2.26 filters, as well as less pronounced 
or non-existing lags near vi s ual mi nimum, confirming 
the findings of ISmith et al.l (|2002), as can be seen in 
Figure [3] The reported phase shift of w 0.15 - 0.22 
with which the NIR maxima lag behi nd visua l max- 
ima typical for M -type Miras (Nad zhip et al.l l200ll 
Smith et alJ2002 ) is also clearly s een in our da t a. Als o 
in accordance with the findings of lSmith et al.l d2002h . 
we find that the maxima in longer wavelengths (i.e. in 
the L band) occur before the K maxima but after the 
visual maxima for all of our objects except for W Hya, 
where the L light curve is too peculiar to be fitted with 
a sine function. 

Both Visible and NIR lightcurves given in Figure [3] 
span the observation period of each object and clearly 
illustrate the substantial differences between the cycle 
amplitudes and pulsation periods. 

3.2. Multi-wavelength, multi-epoch UD angular 
diameters 

Figures H] to [6] show the UD angular diameters of 
the sample of Miras at various visual phases (see Ta- 
bles @] to [9]), folded with phase, together with visual 
photometry data from the AAVSO (A.A. Henden et al. 
2006, private communication). 

Near-infrared pulsations of Mira stars h ave been ob- 
served in terferometric a lly bef or e, e.g. by 



Tuthill 



etal.1 

dl995h: IPerrin et al.l dl999l): lYoung et al.1 (bOOOh: 



Thompson et al.1 d2002l); I Woodruff et alJd2004l) ; lFedele et al 
d2005l)lRagland et al.l d2006h . thoughnever with the 
wavelength, phase and cycle coverage of this study. 
We have detected diameter variations for all Miras in 
this study. As discussed for the light-curves in Section 
3.11 the physical diameter vs. phase curves are cer- 
tainly not strictly sinusoidal. Nonetheless, a sinusoidal 
shape of the curve can be observed, and although not 
perfect, the sine curves fitted to the data have been in- 
cluded to guide the eye. 

With the extensive coverage of phases and cycles in 
four filter bandpasses, it becomes possible to probe the 
stellar atmosphere both for geometric pulsation of the 
continuum forming layers (the so-called photosphere), 
and for contamination of the continuum by molecular 
blanketing. This allows us to further constrain existing 
theoretical models and to make more sophisticated de- 
mands on future models. R Hya will be exempt from 
further discussion as the data sampling is too sparse 
(see Figure|6] lower panel). 



4. Discussion 




1.0 1.2 



Wavelength [^m] 



Fig. 7. — Near-infrared spectrum of Mira R Cha at two 
different pulsation phases (1996Ma y26 at phase 0.6 
and 1 996March3 at phase 0.3) from lLancon & Wood 
(2000), showing the position and shape of the z 1.08, 
J 1 .24, H 1 .65 and K 2.26 filters. The J 1 .24 and H 1 .65 
filters penetrate to layers that lie closer to the contin- 
uum forming photosphere, whereas the K2.26 filter 
sees a portion of the spectrum that is more contami- 
nated by molecular opacities. 



4.1. Effects of Molecular Absorption on Multi- 
wavelength Diameter Observations 

Figure|7]shows the NIR spectrum of a Mira (R Cha) 
at two different pulsation phases, with the position and 
shape of the z 1.08, J 1.24, H 1.65 and K2.26 filters, 
and Figure [8] shows the ISO spectrum of R Cas, with 
the position and shape of the L 3.08 and L 3.31 filters. 
The narrow band filters J 1 .24 and H 1 .65 should show 
little molecular contamination by absorption bands, 
i.e. they are close to sampling the continuum-forming 
layers (see, e.g., 
tamination effect 



Tei et al.l 120031) . whereas the con- 
should be more prominent in the 
K2.2 6 and in particular in the L-band filters (see 
e.g., Jacob & S chol 
Ireland et al .112004 



o3 

a.b) 



2002; Mennesson et al. 2002; 



Within the L band, our filters (L 3 .08 and L 3 . 3 1 ) sam- 
ple portions with different molecular absorption lines, 
origi nating mainly from H? Q (but also from OH and 
SiO. lMennesson et alJ2002l) in Mira atmospheres. The 
L3.08 bandpass lies deeper in the absorption feature 
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Fig. 4. — Uniform disk angular diameters for o Cet (upper panel) and R Leo (lower panel), folded with pulsation 
phase. The symbols with error bars show UD diameters for the J 1 .24, H 1 .65, K 2.26 and L 3.08 filters and the dotted 
lines show the best-fitting sinusoids. The dots show the visual magnitude from the AAVSO database averaged into 
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Fig. 8. — In frared ISO spectrum of R Cas at various 
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ble[2j overlaid. The stronger absorption, mainly due to 
H2O molecules, in the L3.08 filter comes from cooler 
strata (and thus further away from the continuum form- 
ing layers) and leads to the measurement of greater an- 
gular diameters. See text (Section l4~H 



than the L3.31 bandpass. The photons seen through 
the L 3.08 filter should therefore originate from cooler 
strata that lie higher in the stellar atmosphere than the 
L3.31 photons. 

Figure [9] shows the diameter vs. wavelength rela- 
tionship for o Cet at four different phases. It clearly 
shows the effect of varying absorption features on the 
perceived diameter. As our filters coincidentally sam- 
ple more opaque (contaminated) layers with increas- 
ing wavelength, there is a perceived diameter increase 
as a function of wavelength, with the exception of the 
J 1.24 and L3.31 filters, which lie further out of the 
molecular absorption bands sampled by the z 1 .08 and 
L3.08 filters, respectively. This increase of diame- 
ter with wavelength is obvious throughout the pulsa- 
tion cycle for o Cet, all other Miras showing excur- 
sions from the J 1 .24 > H 1 .65 > K 2.26 diameter trend 
at some phases. In particular R Leo seems to show 
a much more complex layering (see Section 1431 with 
Figure|4) . All stars in our sample show a large (> 30%) 
increase in UD diameter between the K2.26 filter and 
the L3.08 fi l ter, firs t noted in the case of R Aqr by 
iTuthill et al.l (I2000al) . This increase indicates that the 
L3.08 filter samples a molecular layer at a consider- 
able distance from the photosphere. The distance be- 
tween the layers varies with phase and will be fur- 
ther discussed for each object individually, as it differs 
greatly from star to star. 
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Fig. 9. — Diameter vs. wavelength relationship for 
o Cet measured at visual phases 0.1, 0.3, 0.7 and 0.9. 
All visual phase measurements carry a binning error of 
±0.1 



4.2. Phase dependence of Multi-wavelength Near- 
Infrared Diameters 

The distinction between true photospheric pulsa- 
tion, i.e. the upward and downward motion of the 
continuum forming layers, and the effects of molec- 
ular blanketing is not trivial to unravel, even with high 
angular resolution data. The varying molecular opac- 
ity is dictated by density changes as a shock front 
travels through the stellar atmosphere and temperature 
changes due to variations in the radiation field, causing 
molecules to dissociate and re- form (cf. IScholj|2003l 
and references therein). These more or less opaque 
strata which vary with the pulsation phase (and be- 
tween cycles) can make it difficult to derive the near- 
UD diameter of the underlying geometrically pulsat- 
ing continuum layer, veiling the geometric amplitude 
and the phase of pulsation. As the brightness distri- 
butions of Miras are more complex than simple UDs, 
changes in molecular opacities of different layers may 
be capable of mimicking changes in our derived diam- 



19 



eters, complicating the interpretation of our data. The 
behaviour of the S-type Mira x Cyg differs from the 
other Miras in our sample in various respects and will 
be discussed separately. 

Of the three filters within the JHK bands, the at- 
mospheric opacity reaches its highest values in the 
contaminated (albeit less contaminated than the stan- 
dard K bandpass) filter K2.26 throughout most of 
the pulsation cycle, as evidenced in Figures [4]|6] by 
the large angular diameters. This increased molec- 
ular o pacity has been predic t ed by various mode ls 



the relative diameter pu lsation amplitude for S Lac of 
Thompson et al.l (120021) . who report a 12%-21% peak- 



fe.g. IJacob & Schold Eool llreland et all l2004allbl 



Ireland & Schold2 006) and has been observed by, e.g. 



Millan-Ga bet et al.l (2005). 



A diagnostic observable from our data which is 
more robust against sources of systematic errors (e.g. 
change in seeing, wind induced wobble of the tele- 
scope, calibrator characteristics) is the relative diame- 
ter variation. Figure[lO]shows the relative UD diameter 
variation as a function of phase for the 5 M-type Miras 
in this study, averaged into phase bins. The relative 
UD diameters were obtained by dividing each object's 
measured angular diameters by its mean, and the er- 
rors are representative of the scatter within each phase 
bin. 

The data illustrate the homogeneity within this sample 
of Miras, and allow us to investigate systematic differ- 
ences between the observed bandpasses. In order to 
extract phase offsets and relative pulsation amplitudes 
in each filter, we fitted sinusoid cycloids to the com- 
bined data, without implying that the objects vary in 
such a simple fashion. 

The ensemble shows a systematic phase-shift of the 
relative UD diameter variation with wavelength, not 
unlike the trend reported for the light-curves in sec- 
tionO The J 1.24, H 1.65, and K 2.26 UD diameters 
reach their minimum at approximately phases 0.9, 1.0, 
and 1.1 respectively. The L3.08 bandpass diameters 
behaves differently, and will be examined in more de- 
tail at the end of this section. This disagrees with th e- 
oretical model predictions by Ireland et al. ( 2004afbh . 
where the diameter minimum should occur at roughly 
phases 0.7-0.8 for all bandpasse s , but i s comparable to 
the findings of IThompson et alj ([2002), who followed 
the M-type Mira S Lac through its pulsation cycle in 
various sub-filters of the K band. Using the sample 
as a single "artificial" Mira also allows us to verify 



to-peak sinusoid amplitude in the Tif-band, with our 
value of 14% peak-to-peak pulsation in the K2.26- 
band. The peak-to-peak sinusoidal pulsation ampli- 
tudes for the J 1.24, H 1.65 and L3.08 bandpasses are 
14%, 22% and 6% respectively, also disagrees with 
theoretica l models which pre dict much higher ampli- 
tudes (cf. llreland et al.ll2004allbh . 

According to current models (e.g. llreland et al. 



2004a b), more contaminated layers (such as the ones 
sampled by the K2.26 and L3.08 bandpasses) should 
experience slightly greater diameter pulsation ampli- 
tudes. 

This trend cannot be seen in our observations: the 
layer experiencing the greatest diameter pulsation is 
the less contaminated H 1 .65 layer. Explanations for 
this could be: 

(i) The relatively narrow width of our K2.26 filter 
compared to the standard K filter used in the models, 
combined with the position of our central wavelength 
in a possibl e minimum of molecula r contamination as 
reported by Thom pson et alJ d2002h . 

(ii) A periodic change of optical depth of a layer, offset 
in phase to the photospheric pulsation, could result in 
large variations of the observed pulsation amplitude. 

(iii) The assumed spherical symmetry could be signifi- 
cantl y violated in the ou ter, more contaminated, layers 
(cf. iRagland et al . 2006 and references therein), effec- 
tively _cauj>mgjle£artu^^ (cf., 
e.g jHofmann et alj|2000lllreland et aU2004cl) . 

(iv) The theoretical models need revision to accommo- 
date for these observations. 

Higher molecular opacity can be expected at near- 
minimum phases, when the outer layers a re cooler and 
more molecules are formed (see, e.g., llreland et al.l 
2004a,b). For our sample stars (with the exception 
of x Cyg), the maximum apparent JHK UD angular 
diameter values are typically found near minimum vi- 
sual phase, supporting existing model interpretations. 

A diagnostic observable from our data which is 
more robust against sources of systematic errors (e.g. 
change in seeing, wind induced wobble of the tele- 
scope, calibrator characteristics) are the relative diam- 
eter ratios. FigureQT]shows the UD diameter ratio be- 
tween the different filters for all our sample stars plot- 
ted vs. phase. There is no obvious dependence on pul- 
sation phase of the H 1 .65/J 1 .24 quotient, though there 
is noteworthy scatter for our sample of 5 Miras (R Hya 
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Fig. 10. — Relative angular diameter variation as a function of phase for the 5 M-type Miras in this study. The sine 
curves represent fits to the ensemble of angular diameters divided by the mean angular diameter for each star. The data 
has been averaged into phase bins and is repeated to show two cycles for better recognition of the pulsation pattern. 
The results for the different filters are offset for clarity. 
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was never observed simultaneously in both filters) that 
might mask a minor phase effect . The mean value of 
this diameter ratio is 



Angular Diameter Ratios 



R 



H 1.65/1 1 



.24 = 1.02 + 0.10 



(1) 



which agrees, to withi n errors, with the value of R q/j = 
1 .08+0.09 reported by lMiilan-Gabet et aD d2005h . The 
lack of a phase dependent signature indicates a close- 
ness in temperature and opacity variations, as can also 
be derived from the closeness in phase and pulsation 
amplitude seen in Figure [10] The diameter ratio close 
to unity shows the geometric closeness of the two lay- 
ers. 

Note that of all stars, o Cet is the only one that has 
a ratio that includes some phase dependent effects and 
is slightly smaller than unity. 

The ratio between the K2.26 and H 1.65 filters has 
a more pronounced pulsation phase signature (see Fig- 
ure QT| center panel), which reflects more complex 
and disjoint temperature and opacity changes between 
these two layers. The ratio reaches its minimum (i.e. 
H 1.65 UD > K2.26 UD) before minimum light and 
its maximum at maximum light. I n order to compare 
our da ta with the observations of iMillan-Gabet et al.l 
(12005b , we also calculated the mean diameter ratio to 
be 



_< s . 




#K2.26/H1.65 = l.H ±0.11, 



(2) 



a mean ratio marginally larger that the H 1.65/J 1.24 
mean ratio. Again, our data agree, to within errors, 
= 1.12 ± 0.09 reported by 



with the value of R k 



Millan-GabetetaLl (2005). 



/H 



As mentioned earlier, the L3.08 angular diameter 
behaves differently. It reaches minimum values at min- 
imum light and the best fit to the diameter pulsation is 
shifted by 0.5 cycle compared with the H 1.65 pulsa- 
tion, albeit at a much smaller relative pulsation am- 
plitude. As the lower panel of Figure [TTl shows, this 
phase shift and the UD diameter ratio between the two 
layers seems to be similar for all Miras in our sample 
(including^ Cyg). For consistency sake, we calculated 
the mean diameter ratio to be 



R 



L3.08/H 1.65 



1.81 ±0.24. 



(3) 



This unusual UD diameter variation has never been 
observed before and raises questions about the mecha- 
nism of the observed pulsation in the L 3.08 layer. The 



1.0 
phase 



o Cet A - 
R T.co O 
li Cas x - 




Fig. 1 1 . — Angular diameter ratios for all Miras in this 
study. The solid lines indicate the mean ratio and the 
dashed lines indicate the standard deviation of the sam- 
ple. The data is repeated to show two cycles for better 
recognition of the pulsation pattern. 
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L3.08 light-curve follows a similar trend to the JHK 
light-curves, possibly indicating a temperature devo- 
lution similar to the lower layers, and as the inner lay- 
ers of the star are shrinking and heating up, the outer 
layers are either expanding or becoming increasingly 
opaque. Whether opacity effects or dynamic motion 
of these outer layers (or both) are responsible for this 
surprising behaviour is the subject for model interpre- 
tations and will be the subject of a subsequent study. 

4.3. Individual Stars 

Of the 6 Miras observed, 4 have observational phase 
coverage suitable for further comparisons with pulsa- 
tion models (o Cet, R Leo, R Cas and x Cyg). The 
Miras studied differ substantially in behaviour regard- 
ing pulsation amplitudes, diameter-wavelength rela- 
tionships and diameter-phase relationships 

In this section we discuss the results of this study 
for each individual star, emphasising the main differ- 
ences and similarities found in this subset. 

o Cet 

The prototype of Mira stars is one of the most ob- 
served variable stars, due to its brightness, amplitude 
(Vx! 10 - 3, A. A. Henden et al. 20 06, private commu - 
nication) and closeness ( 1 07 +6 pc, iKnapp et alJ2003h 



Its size in diff erent bandwidths (e.g. 



Haniffet al. 



1995 jMennesson et all l2002MWoodruff et alJ 12004). 



optical spectr a Jjovl Il954 j ). lightcur v es in different 
color s (e.g. IWhitelock et alJ uOOOt Nadzhip et al 
2001; AAVSO), a symmetries (e.g. iKarovska et al 
1 00 " iTuthill et al.1 Il999[). and companion star (e.g. 



2001 
1991 

Karo vska et alJl997UWood & Karovskal2006llrreland et al 
2007) have been subject to intense research. 



Figure H] shows the UD diameter variation in the 
J 1.24, H1.65 and K2.26 bandpasses as nearly syn- 
chronous (within the 0. 1 phase shift shown in Figure 
[Tol l, sinusoidal pulsations, with an apparent phase shift 
to the diameter pulsation in the L 3.08 bandpass. The 
UD diameter vs. phase curves agree well with the fit- 
ted sine functions, with reduced^ 2 for J 1.24, H 1.65, 
K2.26 and L3.08 having the values 1.45, 1.17, 1.02 
and 0.80 respectively. The shift of sa 0.5 cycles be- 
tween the sine curves fitted to the diameter vs. phase 
in the H 1 .65 and the L 3.08 bandpasses has never been 
observed nor predicted, and can also be observed in R 
Leo, R Cas and W Hya. A similar shift, though not as 
pronounced, can be seen in^ Cyg 



The UD diameters vary between 22.2 + 1.0 and 
29.9 ± 1.9 mas (« 35%) in the J 1.24 bandpass, and 
between 53.2 ± 1.4 and 61.5 + 2.3 mas (a 16%) in 
the L3.08 bandbass. The K2.26 UD angular diam- 
eter values lie in the range of 31.0 + 1.9 mas near 
V- maximum and 37.0 + 1.2 mas at phase 0.7. These 
values agree with the interferometric UD diameters 
within th e K-band of | Ridgwav et alJ (119921) (phase 
0.8) and IWoodruff et all d2004h (phases 0.1-0.4, see 
Table [TDl . Our UD angular diam eters are gener- 



ally la rger than those measured by iMennesson et al. 



d2002h (phases 0.9-0.0). This is most likely due to to 
the larger spatial frequencies (i.e. longer baselines) 
at which their measurements where made, combined 
with the known dep arture of o Cet's brightness distri- 
bution from UD (cf. IWoodraff et a"i1l2004l) . 
Our K2.26 UD diameters are also large r than the 



molec ular layer diameters obtained by iPerrin et al. 
d2004h (24.95+0.10 to 26.84 + 0.06 mas) at similar 
phases by fitting ad hoc scenarios (a photosphere sur- 
rounded by an emissive and absorbing layer) to K and 
L' interferometric data. This is to be expected, as their 
visibilities also show obvious departures from simple 
models, and fitting the same data to a brightness dis- 
tribution consisting of a central object with a bright 
molecular shell would yield smaller diameters than a 
UD fit. Note that these very simple scenarios are not 
always unique and a new ad hoc parameter set has to 
be determined for each observation, making a com- 
parison awkward, o Cet shows little cycle-to-cycle 
variation throughout our data, although this could be 
attributed to the observing of a stable era of its pulsa- 
tion, and might have been different if observed, e.g., 
10 years earlier. 

The L3.08 » K2.26 > H 1.65 > J 1.24 layering 
of monochrome diameters is strictly monotonous, un- 
like, e.g., R Leo, where the sequence of UD diameters 
seems to invert during the pulsation cycle. The fact 
that o Cet appears largest in the JHK filters approx- 
imately at visual phase 0.6, roughly coinciding with 
the minimum of the NIR and visual light curves (see 
Section [3~TT i is in accordance with various model in- 
terpret a tions (e.g. Ijacob & Scholzl 12002 llreland etal 



2004a b: Irelan d & Scholzl l2006h. The unusual be 
haviour of the L3.08 layer, as described in Section 
14.21 will be further discussed in Section [5] and in a 
follow on paper. 



23 



Table 10 

Published interferometric UD angular diameters in the K band for the stars in this study 



olar 


Reference 


Angular 


Centre Wavelength 


.Bandwidth 


Visual 






I |1 n m at at* I tTI Q c I 

LVldlllcLcr [IlldSJ 


[^m] 


[/zm] 


l-Jr* OCA 

r Ild.Sc 


o Cet 


Ridgwav etal. (1992) 


36.1 + 1.4 


2.2 


0.4 


0.8 




Mennesson et al. (2002) 


28.79+0. 10 


2.2 


0.44 


0.9 






25.73+0.09 


2.03 


^0.1 


0.0 






25.13+0.08 


2.15 


^0.1 


0.0 






25.19+0.12 


2.22 


* 0.1 


0.0 






29.22+0.12 


2.39 


* 0.1 


0.0 






24.40+0. II 


2.16 


0.32 


0.0 




Woodruff et al. (2004) 


29.24+0.30 


2.2 


0.4 


0.1 






29.53+0.30 


2.2 


0.4 


0.2 






30.49+0.30 


2.2 


0.4 


0.3 






33.27±0.33 


2.2 


0.4 


0.4 


R Leo 


di Giacomo et al. (1991) * 


33.0+1.3 


2.16 


0.03 


0.2 




Perrin etal. (1999) 


28.18+0.05 


2.16 


0.32 


0.2 






30.68+0.05 


2.16 


0.32 


0.3 




Tei et al. (1999)* 


34±2 


2.36 


0.46 


0.4 




Monnier et al. (2004) 


30.3+0.3 


2.26&2.16 


0.05&0.32 


0.6-0.7 




Millan-Gabet et al. (2005) 


29.91 +0.27 


2.16 


0.32 


0.4 




Fedele et al. (2005) 


28.1+0.05 


2.2 


0.4 


0.1 






26.2+0.01 


2.2 


0.4 


0.0 


RCas 


van Belle et al. (2002) 


22.03+, 2 , 1 / 


2.2 


0.4 


0.8 


WHya 


Monnier et al. (2004) 


42.5+0.7 


2.26&2.16 


0.05&0.32 


0.5 




Millan-Gabet et al. (2005) 


39.9+0.2 


2.16 


0.32 


0.6 


*Cyg 


Mennesson et al. (2002) 


23.24+0.08 


2.16 


0.32 


0.38 


RHya 


Millan-Gabet et al. (2005) 


23.9+0.5 


2.16 


0.32 


0.8 






25.8+0.2 


2.16 


0.32 


0.8 



'obtained by lunar occultation 
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RLeo 



least squares sinusoidal fit (reduced^- 2 = 1.38). 



Although R Leo has a similar period (310 days), V 



magnitude range (4.4 - 11.3, iKholopov et alj 11998) 



and spectral type as o Cet,we found significant differ- 
ences in this study. The pulsation amplitudes of layers 
in different bandpasses and their phases with respect 
to one another seem to show a more complex trend 
than o Cet's atmosphere. The pulsation is less pro- 
nounced than o Cet's in the NIR, with UD diameters 
varying between 29.6 ± 1.4 and 33.1 +2.0 mas (« 12%) 
in the J 1 .24 bandpass, which shows the largest rela- 
tive pulsation amplitude. The cycle-to-cycle variation 
in the diameter pulsation is significant in the K2.26 
bandpass (compared to, e.g., o Cet). Again, this effect 
could be due to mere coincidence, and only long term 
(> 20 years) observation campaigns can ascertain this 
variability. This makes it difficult to compare our K 
UD diameter measurements with values found in the 
literature (see Table [TOt. There is nonetheless an over- 
all agreement with cited UD radii except for the lon g 
baseline UD angular diameters of iFedele et al.l (|2005). 
This can be due to the same effect as described in Sec- 
ti qn o Cet concerning the l ong baseline measurements 



of Mennesson et al 



( 20021) . With our H 1.65 measure- 
ments we are able to disambiguate the two possible 
solutions (due to ambiguo us model fitting) fo r the H - 



band diameter found by iMillan-Gabet et al.l ([2005) 



Due to the seeming lack of cycle-to-cycle dependence 
in this bandpass, we can reject their smaller angular 
diameter of 23.8±0.3mas at phase 0.4 and substanti- 
ate their larger diameter of 32.4+0.4 mas 

The variation between 47.8 ±1.8 and 56.5 ±2.0 mas 
(a; 18%) in the L3.08 bandpass shows a steady in- 
crease in diameter between 97Decl6 and 01 Junl 1 and 
no significant diameter change up to 3 pulsation cycles 
later. This gradual increase in angular diameter over 
4 cycles may be understandable in terms of the non- 
cyclic time evolution of positions of outer mass zones 
of pulsation models over several successive cycles (cf. 
Figures 1 of Ireland et al. 2004a, b), which affect the 
position and physics of water shells. This interpre- 
tation is supported by the steady brightening of the 
visual photometry maxima during the same 4 cycles 
in which the UD diameter increases, indicating vari- 
ability in timescales longer than the pulsation phase. 
The pulsation signature of the L 3.08 bandpass shows 
nonetheless the general w 0.5 phase lag compared to 
the H 1.65 band, albeit with a larger uncertainty in the 



R Cas 

Of the two Miras with the longest periods in our sam- 
ple, R Cas (430 days) also exhibits the latest spectral 
type. In fact, at minimum visua l pulsation phase it 
define s the spectral type M10 (see Lockwood & W ing 
197 II) . It seems to follow the model-predicted phase- 
diameter trend of larger diameters around minimum 
light in the H 1 .65 filter, albeit with some cycle-to- 
cycle uncertainty. The scatter is considerably greater 
in the J 1.24 and K2.26 filters, but the same trend is 
still present. Note that the relative phase shift of di- 
ameter maxima from shorter to longer wavelengths in 
JHK is very distinctive, even when the uncertainties 
introduced by cycle-to-cycle variations in the diame- 
ter pulsations are taken into account. The position of 
the L3.08 layer also follows the pattern of the other 
M-type Miras, in that it seem s to pulsate with a 0.5 
phase shift to the H 1.65 layer. Ivan Belle et alj (2002) 
measured a K UD angular diameter of 22.03* 2 ," mas at 
phase 0.81, which does not agree with our diameter of 
28.9±1.9 at the same pulsation phase but agrees with 
the diameter measured at roughly the same phase in a 
different cycle (26.1 ±2.1 at phase 0.74), another indi- 
cator of noticeable cycle-to-cycle variation. 



WHya 



At a distance of 78±3 pc dKnapp et alJl2003l) . W Hya 
is the closet and best resolved Mira in our sample. 
We observe a relatively time independent K2.26 an- 
gular diameter of * 40 ± 5 mas in the phase range 
O = 0.53 - 0.88. Because of W Hya's low declination 
and proximity to the ecliptic, and because the period 
is close to 1 year, the phase coverage only spans half 
the pulsation cycle (see Table |7|. In this short time 
interval, covering only 0.44 phases, the J 1.24 diam- 
eter shrinks from 47.1 ±2.7 mas to 32.0 ± 2.3 mas, a 
32% decrease, whereas the H 1.65 diameter decreases 
by 22% and the K2.26 diameter shows a decrease of 
*20%. 

The previ ously measured UD an gul ar diameters in the 
K-bm d by Mo nnier et"ai] d2004 and lMillan-Gabet et al. 
d2005l) (see Table ITQb are consisten t with our measure- 
ments. The //-band diameter of Millan-Gabet et al. 



(|2005l) (31.3±0.3mas) is slightly smaller than our de- 
rived H1.65 UD angular diameters, which is most 
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likely due to the same effects of long baseline interfer- 
ometry as described for o Cet and R Leo. Although 
all Mira intensity distributions show deviations from a 
UD profile, we can clearly detect those only in W Hya. 
This could be due to the presence of dust emission 
and/or the partially resolved brightness distribution of 
extended molecular layers in the upper atmosphere. 

The largest L 3.08 angular diameter is found around 
maximum visual light, where the J 1 .24 and H 1 .65 di- 
ameters are smallest. 

X Cyg 

The Mira in our sample with the second longest 
pulsation period, x Cyg (408 days) is the only S- 
type Mira ,albeit with M-type characteristics (cf. 
Keenan & Boeshaarll980l) . It deviates noticeably from 



maximum around minimum light, as observed in the 
other Miras in our sample, suggesting a deeper view 
into;^ Cyg's atmosphere. 

Young et al. (2000) only detected a slight variation 



the M-type Miras in many ways. The visual light- 
curve of x Cyg is well known for its particularly large 
magnitude range, and its nea r-IR colours are also al l 
unusually red, as observed by IWhitelock et a l. (2000). 
The L3.08 mean magnitude, in particular, is compa- 
rable to or brighter than that of K2.26, whereas for 
all other Miras in our study the L 3.08 magnitudes lie 
in the rage between the H 1.65 and the K2.26 light- 
curves. 

The relative diameter pulsation amplitude is larger 
in all filters than that of the other Miras in this pa- 
per, and the derived UD angular diameters show small 
cycle-to-cycle variation. There are no model predic- 
tions for S-type Miras to date, and we can only spec- 
ulate whether or not this behaviour is related to the 
fundamental stellar parameters or to the fact that dif- 
fering C to O ratio lead to si gnificant chang es in H2O 
formation characteristics (cf. lOhna ka 2004) and more 
stable opacity structures. 

While the L 3.08 UD diameter pulsation follows the 
same trend as o Cet, R Leo, R Cas and W Hya, in that it 
reaches maximum amplitude around maximum visual 
light, the JHK diameters show a somewhat different 
behaviour. Note that the shapes of the visibility curves 
deviate increasingly from a UD profile with increasing 
wavelength. The K2.26 and H 1.65 UD angular diam- 
eters reach their maximum values just after maximum 
light, at phase 0.1-0.2, coinciding with the maximum 
of the correspondent light curves. This may indicate a 
heavy contamination of these filters by molecular lay- 
ers in the star's atmosphere. On the other hand, the 
J 1 .24 band data infers a smaller UD diameter with its 



in the /-band angular diameters, albeit over a very 
small phase coverage. If we fit a sine curve to our 
derived J 1 .24 angu lar diameters, then we find that 
Young et all j2000l) measurements lie roughly sym- 
metrically around the diameter pulsation minimum, a 
factor that could explain the lack of diameter variation 
in their studie s. We converted the /-band Gaussian 
FWHM from lYoung et alj d2000b to UD angular diam- 
eters usin g the UD to Gau ssian ratio R « 1.5 calcu- 
lated by iBurns et all dl998l) . and find that the diame- 
ters appro ximately match our sim ple sine curve pre- 
dictions. Me nnesson et alj d2002l) report a K' (cen- 
tered at 2.16^im, 0.32yum wide) UD angular diameter 
of 23.24±0.08mas for^ Cyg on 2000 May, at variable 
phase 0.38. On June of the same year (phase 0.48) 
we obtained a K2.26 diameter of 25.0+2.2, in reason- 
able agreement. In the next two subsequent years (both 
at phase 0.36) we observed UD diameters of approx- 
imately 30 mas, an increase in diameter of 20% over 
one cycle, which might be linked to transient opacity 
structures. 

5. Summary 

We have measured the diameters and NIR light- 
curves of 6 Miras stars at up to 19 separate phases in 4 
filters, the first study of this magnitude. We present the 
first narrowband 3.08 /vm light-curves of Miras. The 
NIR light-curves can be approximated by a sine func- 
tion and confirm the phase shift of » 0.15 - 0.22 by 
which the NIR maxima la g behind the visual max- 
ima previously rep orted by iNadzhip et alJ (1200 II) and 
Smith et alj d2002l) . In addition we find a NIR photo- 
metric flux relation: J 1 .24 < H 1 .65 < K2.26, in agree 



ment with the JHK observations of IWhitelock et al. 



d200d) and ISmifh etaD d2002l) . The L3.08 fluxes in 
our sample are slightly less or equal to the K2.26 
fluxes, except for^ Cyg, where the L3.08 magnitudes 
are of the comparable to the ones in the K2.26 band- 
pass. 

We found no correlation between NIR photometry 
and UD diameter cycle-to-cycle variations, yet found 
some correlation between the L 3.08 UD diameters and 
the visual light-curves (see section R Leo). 
All observed stars show variations of their UD angular 
diameters as a function of pulsation phase. We find 
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the UD diameter relation J 1.24 < H 1.65 < K2.26 to 
be an average value only, with deviations throughout 
the pulsation cycle, revealing the complexity of phase 
dependant opacity contamination from molecules in 
different layers. Of the Miras in our sample, only 
o Cet shows this strict layering throughout its cycle, in 
agreement wit h theoretical models d esigned to repre- 
sent o Cet (see I Jacob & Scholzll2002h . 
The pulsation amplitude also does not follow the 
model-predicted dep endency on molecular opacity 



([Ireland et al. l2004alTbh . The layers exhibiting the 



largest relative variation in UD diameter (6%-18%) 
were those seen through the H 1.65 and J 1.24 filters, 
which should display less molecular contamination 
and thus less diameter variation. The K2.26 layer has 
UD diameters pulsation amplitudes between 4% and 
7%. 

The molecular layer probed by the L3.08 bandpass 
is significantly further from the photosphere than the 
JHK layers. It varies in UD angular diameter by as lit- 
tle a 4% and as much as 8%, with a 0.5 phase offset to 
the H 1.65 pulsation, and is between 1.5 and 2.5 times 
larger than the H 1 .65 UD angular diameter, depending 
on the pulsation phase. This behaviour has not been 
observed before, and detailed models are needed to 
understand it. 

Another quantity predicted by models are the di- 
ameters in different wavelengths at different phases. 
When comparing the two stars whose parameters are 
thought to resemble the model input parameters the 
most, R Leo a nd o Cet (both wi th revised Hipparcos 
parall axes, see KnarJD et al |2003l) . with model predic- 



tions (llreland et al. I l2004allbh 7 we find that our UD di- 



ameters are significantly too large in the K2.26 band- 
pass. This could be explained by too high model ef- 
fective temperatures, but a more careful, model based 
interpretation is needed to understand this effect. 
Given the small baselines (< 10 m) used in this exper- 
iment, it is remarkable how consistent the data is for 
the M-type Miras. It is even more surprising how dif- 
ferent the S-type Mira^ Cyg appears to be when it's 
mutli-wavelength pulsation signature is analyzed. The 
H 1.65 and K2.26 UD angular diameters are small- 
est around minimum light, contrary to all models and 
previous observations of M-type Miras. Further work 
with these data, including model comparisons with in- 
dividual stars, imaging and asymmetry studies, will be 
presented in subsequent publications. 
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